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1. Introduction 

There are a number of issues and questions about baryons that can be addressed using lattice 
QCD. For example the level ordering between the mass of the parity partner of the nucleon and 
the mass of the first excited state of the nucleon. Also as we discuss in section || there are matrix 
elements of the nucleon, such as the scalar current of strange quarks in the nucleon that have an 
uncertainty of at least a factor of 3. In this paper we make a start to study the mass of the parity 
partner of the nucleon and develop lattice techniques to compute the strangeness content of the 
nucleon using staggered fermions. 

Seduced by the light sea quark masses and large volumes of the unquenched lattice QCD 
calculations performed by the MILC collaboration [|l|, 0, we use the improved staggered fermion 
formulation (known as "Asqtad") to do our calculations. The unquenched calculations use 2 + 1 
flavours of sea quarks with quark masses of 0.007/0.05. The lattice volume was 20^ 64 and P = 
6.76. We used the gauge configurations generated by the MILC collaboration stored on the public 
NERSC archive. 

2. The effect of open decays on the lattice. 

Many of the really interesting questions in hadronic physics, such as the existence of hybrid 
mesons or glueballs, requires lattice QCD to be able to cope with hadrons that decay via the strong 
interaction. Recent unquenched lattice QCD calculations now have light enough sea quarks that 
hadrons are starting to decay. The problem is most severe for particles that decay via S-wave de- 
cays, because particles that decay via P-wave channels have additional kinematic obstacles because 
the decay products have to have momentum and this can be quite large on typical lattices 

Indeed the MILC collaboration have claimed to see problems for the state [|l|, |2|] and 1 ^ 
hybrid meson [Q], due to the quark masses being light enough for these particles to decay into two 
mesons. There have also been problems reported with the singlet 0++ mesons [||]. The experimen- 
tal status of the above mesons is confusing, so it would be better to study the effect of hadronic 
decays on hadrons where there is less uncertainty about the experimental data. So it is important to 
study the issue of open decay channels of baryons. 

The masses of baryons quoted in the particle data table (PDG) are usually extracted using a 
partial waves analysis, hence the angular momentum of the decay can be read off from the quantum 
numbers. In the table ^ we summarise some data about some baryons from the PDG [^. From 
table ^the decay of the parity partner of the nucleon (N(1535)) looks like a good candidate to have 
the open decay of Nn, in the current generation of unquenched calculations. 

There has been a lot of work on using quenched QCD to study excited baryons. It has been 
found that large volumes (> 3 fm ) are required, as well as light quark masses for the chiral extrap- 
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olation [^. At the moment only the improved staggered program of the MILC collaboration comes 
close to these requirements. Although in this calculation we have used a lattice with a physical box 
side of 2.4 fm, the MILC collaboration have also generated unquenched gauge configurations 
with a box size of 3.4 fm. 

Unfortunately, the identification of baryon operators with physical states is non-trivial in the 
staggered formalism [^]. The staggered baryon operators that are local in time break down into 
irreducible representations like: 5.8 + 2.8' + 4.16 The nucleon couples to 8 and 16 representa- 
tions [^]. Things are more complicated for the parity partner states. The A(1405)^ will couple to 
the odd 8 representation. There are no interpolating operators for staggered fermions that couple 
to the parity partner of the nucleon in the ground state. 

The staggered parity partner to the nucleon with mass wi_ is extracted from the correlator c{t) 
using the fit model in equation |2.li 

c{t) = a+exp(— m+f) + (— l)'a_exp(— (2.1) 

where is the mass of the nucleon. There have been earlier attempts to extract the physical m_ 
with unimproved staggered fermions [^. 

In figure I] we plot the masses m+ and m_ in lattice units as a function of the starting time in the 
fit. We used the gauge configurations from the MILC collaboration mentioned in the introduction 
and used wall sources in Coulomb gauge to generate the quark propagators. The actual number 
of configurations is in table |[ We also include the published number for m/v_|_ mass from MILC 
(they don't publish numbers for m/v_) as well as the sum of the nucleon mass and the mass of the 
Goldstone pion. 

The data in the figure |I] aix consistent with m/v_ = m^j^ + m,!. However, the box side has a 
physical length of 2.4 fm, so the m^v- niass may be suppressed by finite volume effects. In physical 
units the mass of the Goldstone pion is 310 MeV in this calculation. The physical mass splitting 
between the nucleon and N(1535) state is about 600 MeV, so the state with the opposite parity in 
the nucleon channel is allowed energetically (at least) to decay to a nucleon and pion. A careful 
study of the mp^^ mass in larger volumes and with heavier quark masses (where the decay is closed) 
is required to confirm that this state has decayed. 



We are not very comfortable with using equation 2.1 to extract physical numbers for nij^^, 
because multiple exponential fits to a single channel can be unstable. Also to study the issue of 
particle decay we need as many excited states as possible. So we have started to investigate using a 
basis of interpolating operators (<!>,■(?)) and extracting the masses and amplitudes via a variational 
technique. 

A matrix of correlators is measured 

c{t)^j = {^^{t)^j{Of) (2.2) 
and analysed with a fit model of the form 

c{t)=A^TA (2.3) 

For staggered fermions, the diagonal matrix T contains exponentials {e^'"+') and exponentials with 
prefactors of (—1)'. We prefer to fit the data to the expression in equation rather than use 
techniques such as the generalised eigenvalue method. 



042/3 



Baryon spectroscopy with improved staggered quarks. 



C. McNeile 



Mass of nucleoli and opposite parity state 



1.15 - 

1.1 - 

1.05 - 

1 - 

0.95 - 

= 0-9 - 

0.85 - 

0.8 - 

0.75 - 

0.7 - 

0.65 - 

06 

1 



o 


Mass N+ 


□ 


Mass N-- 


<> 


Mass N+, MILC 




m„+m . 




N pi 



Figure 1: Mass of the nucleon and the opposite parity state that couples to the same correlator as a function 



of the first timeslice in the fit. 2.1 
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Table 1 : Fit results for the nucleon and the state with opposite parity in the same channel. The number before 
the X symbol is the number of time sources for the quark propagator calculations on each configuration 



We used "fuzzed" and local basis states [ ]10| ] originally developed for clover fermions. To 
respect the staggered symmetries fuzzed gauge links only separate the quark and anti-quark by 
multiples of 2. Here we try it for the local staggered nucleon operator that belongs to the 8 repre- 
sentation. 



0{t) 



X even 



(2.4) 



The Chroma software system was used for the lattice calculations [|llj]. 

In table |T] we show some results for the variational analysis versus using wall sources from 
MILC's analysis and our analysis. The variational analysis is inconsistent at the 4a level with the 
results that use single wall sources. For Wilson fermions, wall sources have small statistical eiTors, 
but tend to approach a plateau very slowly. For this data set MILC fitted the nucleon operator at 
smaller times than for the equivalent analysis at other P values. It is not clear what causes the 
difference between the mass from the wall source and variational analysis. 



3. Strangeness content the nucleon 

The matrix element (A'^ | 5.5 1 A'^) is colloquially known as the strangeness content of the nucleon. 
This matrix element is important for baryon chiral perturbation theory. This matrix element is also 
used in the neutralino-nucleon scalar cross-section rtT2l] and hence crucial for the direct detection 
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~0.6 


UKQCD 119|] 


2 
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Table 2: Some lattice and non-lattice determinations of y. 



of dark matter. It is convenient to discuss the ratios of matrix elements 

y = 2 <^l"i^> (3^1) 

{N\uu + dd\ N) 

where A'^ is the nucleon operator. 

In the table ^ we record some recent values for y from lattice and non-lattice methods. There 
has also been a large quenched study of matrix elements related to y by Lewis et al. [|l^]. The value 
of y can range from 0.2 to 0.6. The (A'^ | 55 | A'^) matrix element should not be effected by open 
decay channels, so a precision value should be attainable. 

On the lattice the matrix element (A'^ | | A'^) is extracted from the correlation of a strange loop 
with the two point nucleon correlator [ jn] , p^ ]. The loops were computed via a standard stochastic 
technique using random noise sources [pO|]. With staggered fermions the scalar loop has no phase 
factors [^]. However, the nucleon two point con^elator does have a contribution that oscillates in 
time (see equation ^). We used the unquenched data set described in section |I| with 330 gauge 
configurations. Quark propagators inverted from wall sources were used. 

The strange loop is inserted at time t\ . The nucleon interpolating operators are at times and 



t with t = t\+t2. The fit model for the three point correlator is equation 3.2 



Cf'\h,t2) = £cj'')e-^''';c,-,e-^^'^cf (3.2) 

where a and b label the interpolating operators and the time oscillating term is suppressed. The 
amplitudes and masses M,- are fixed from the fits to the two point function. The (A'^ | 55 | A'^) 
matrix element is proportional to xqq. In these preliminary fits we assumed xq\ =0, but varied xqq 
and x\i. Using fi = 4 and fitting from t=9, we get xqq = 1.5 ± 1.0. For this data set, the two point 
nucleon correlator gets very noisy for ? > 12. Given the large statistical errors we do not quote a 
physical number for (A'^ | 55' | A'^) or y. 



4. Conclusions 



We have started to study baryons on the lattice when the decay channels are open. It is critical 
to include finite volume studies of these states [^. This is computationally feasible with the current 
generation of improved staggered fermion calculations. 

For poorly known quantities such as the strangeness content of the nucleon, we should be 
able to (finally) pin this matrix element down. The statistical error on the bare matrix element is 
disappointingly large, but we hope that by incorporating the variational analysis and using many 
more gauge configurations the errors should be reduced. 
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